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ABSTRACT 

Context. The [Ne II] line 12.81 pm was proposed to be a good tracer of gas in the environments of proto-planetary disks; its origin 
is explained by different mechanisms: jets in outflows, photo-evaporative disk winds driven by stellar X-rays/EUV or by the X- 
ray irradiated proto-planetary disk atmosphere. Previous Spitzer studies gave hints toward the neon emitting mechanism by exploring 
correlations between the line luminosity and properties of the star-disk system. These studies concluded that the origin of the emission 
is likely related to accretion and outflows, with some influence from X-rays. 

Aims. We provide direct constraints on the origin of the [Ne II] emission using high-spatial and spectral resolution observations that 
allow us to study the kinematics of the emitting gas. In addition we compare the [Ne II] line with optical forbidden lines. 
Methods. We obtained high-resolution ground-based observations with VISIR-VLT for 15 stars and UVES-VLT for three of them. 
The stars were chosen for having bright neon emission lines detected with SpitzerflRS. The velocity shifts and profiles are used to 
disentangle the different emitting mechanisms producing the [Ne II] line. A comparison between results from this study and previous 
high-resolution studies is also presented. 

Results. The [Ne II] line was detected in 7 stars, among them the first confirmed detection of [Ne II] in a Herbig Be star, V892 Tau. 
In four cases, the large blueshifted lines indicate an origin in a jet. In two stars, the small shifts and asymmetric profiles indicate an 
origin in a photo-evaporative wind. CoKu Tau 1, seen close to edge-on, shows a spatially unresolved line centered at the stellar rest 
velocity, although cross-dispersion centroids move within 10 AU from one side of the star to the other as a function of wavelength. 
The line profile is symmetric with wings extending up to ~ ±80 km s -1 . The origin of the [Ne II] line is unclear and could either be 
due to the bipolar jet or to the disk. For the stars with VLT-UVES observations, in several cases, the optical forbidden line profiles 
and shifts are very similar to the profile of the [Ne II] line, suggesting that the lines are emitted in the same region. A general trend 
observed with VISIR is a lower line flux when compared with the fluxes obtained with Spitzer. We found no correlation between 
the line full-width at half maximum and the line peak velocity. The [Ne II] line remains undetected in a large part of the sample, an 
indication that the emission detected with Spitzer in those stars is likely extended. 

Key words. ISM: jets and outflows - Infrared: stars - Protoplanetary disks - Stars: formation - Stars: pre-main sequence 

1 . Introduction with a planetary system, we need to better understand the disk 

itself: its composition, size, dynamics, timescales. 



Protoplanetary disks are the outcome of the star formation 
H process. They provide the material for the young star and 

its planetary system. The central star indeed accretes mat- The Spitzer Space Telescope dWerner et all I2004I) made 

ter from the disk until reaching its final mass, after £ a large contribution to the studies of protoplanetary disks 

10 Myr (e.g.. iMannings & Sargent) 119971: lHaisch et ail 1200 it in the mid-infrared (mid-IR) by revealing th e presence of 

iHillen brand 2008; Fedele et al.l20io|) . In order to construct a de- many emission li nes from orga n ic mo le cules (ICarr & N aiita 



tailed picture of the process that turns a pre-main-sequence star | 2008l). water (ICarr & Naiital 120081: ISalvk et al.M 2008 
surrounded by an optically thick disk into a main-sequence star pontoppi dan et all 1201 0|), and atomic species dPascucci et al 

| 2007[ lLahuis et al.l l2007t iFlaccomio et all 12009; Gud el et al 
l2010t iBaldovin-Saavedra et al.ll201 ll) . Among the many emis- 

* Based on observations made with ESO Telescopes Kueyen/UT2 sion lines detected b Y f Ne "] 12 - 81 A«n gained particular 

and Melipal/UT3 at the Paranal Observatory under programmes ID interest. It was proposed to be a good diagnostic of gas in the 

083.C-0471, 084.C-1062, 086.C-0911, and 286.C-5038 u PP er layer of the disk, and is used to study the interaction 

** Now at the CRNS-INSU/UJF-Grenoble 1, Institut de Planetologie between stellar high-energy irradiation and the disk. Currently 

et d'Astrophysique de Grenoble (IPAG) UMR 5274, Grenoble, F- there are three main mechanisms that could explain the presence 

38041, France of [Ne II] in the environment of young stars: 
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- Irradiated disk atmospheres (Glass gold et al.ll2007l) . In this 
scenario, the X-rays from the central star create a warm at- 
mosphere composed of gas in atomic form on top of a cooler 
molecular layer. [Ne II] would be emitted in a region within 
20 AU from the central star, the giant planet formation re- 
gion. 



- Photoevaporative disk wind s driven b y EUV dClarke et al.l 
200UlAlexanderetal.ll2006h. X-ravs (j Ercolano et al.ll2009t 



Owen et al.ll201d) or FUV (iG orti et alJ l2009t). 
- Shock s (IHollenbach & McKed 1 19891 : IHollenbach & Gorti 
2009) and jets dShang et al.ll2010l) . High velocity outflows 
from the star interacting with the surrounding material 
can create strong shocks heating gas to high temperatures. 
[Ne II] would be a good tracer of material at velocities 
higher than 40-50 km s _I . 

Given the modest spectral resolution of the Spitzer infrared 
spectrograph (IRS; R ~ 600, i. e., a velocity resolution of ~ 
500 km s _1 ; iHouck et"aTI 120041) the lines detected are unre- 
solved. Ground-based observations at high spectral resolution 
are needed in order to determine the origin of the neon-emitting 
mecha nism, but the numbe r of high-resolution spectra is still 
small. Hercze g~et al.l d2007l) detected the [Ne II] emission line 
in one out of three spectra of young circumstellar disks ob- 
served with the mid-infrared spectrograph on Gemini North 
(R ~ 30000). Based on the measured linewidth they ruled out 
an accretion flow origin , favoring the photoevaporative theory. 
Ivan Boekel et al.l ([2009) presented high resolution ground-based 
spectroscopy (VISIR-VLT) of the young T Tau triplet, succeed- 
ing in spatially separating the N-S components. These observa- 
tions showed that the [Ne II] emission is strongly domin ated 
by outflows heated by shocks. iPascucci & Ster zik (200§ ob- 
served 6 targets with VISIR-VLT. For 3 stars in the sample, the 
line is spectrally resolved, and the line profiles are consistent 
with those predicted by photoevaporative flow driven by EUV 
from the central star. Pascucci et al .1 (120111) performed an exten- 
sive study on the neon emission in TW Hya using VISIR-VLT. 
The study confirmed the photoevapora t ion as emitting mecha- 
nism of neon in this star. ISacco et al .1 (120121) also studied the 
origin of [Ne II] emission with VISIR-VLT. They observed a 
large sample of young stars and detected the line in 12 out 
of 32 objects. They concluded that [Ne II] emission originates 
mainly in shoc ks for Class I proto stars . This supports the statis- 
tical st udies of iGtidel et al.l (1201 Ol) and Bal dovin-Saavedra et al.l 
(201 I D that jets, i f pres ent, tend to dominate the [Ne II] emis- 
sion. Sacco et al ] (120121) also argued that the emission line stems 
from the inner disk (< 20 - 40 AU) for stars with transition 
and pre-transition disks. Their detailed analysis of the line pro- 
files indicated an origin from a disk wind, although irradiation by 
EUV/X-rays underestimates the blueshift of the line. Our study 
complements the above studies by providing additional VISIR 
spectra for 15 targets, mainly Class II stars, and further com- 
paring the [Ne II] line with optical lines, in particular forbidden 
lines such as [O I], [N II], and [S II]. 



1.1. Sample 

We recently studied the gas emission with Spit z er-IRS 
in a large sample of y oung stars (GudeletaTJ l2010t 
B aldovin-Saavedr a et al. I l20l!t) obtaining a large number of de- 
tections of [Ne II]. For our dedicated VLT- VISIR high-spectral 
resolution study, we selected objects accessible by the VLT 
for which the line fluxes obtained with Spitzer are higher than 
10~ 15 erg cirT 2 s _1 , enough for follow-up ground based ob- 



servations. The stars in our sample are mainly optically thick 
disks (Class II). We included a Class I object ( L1551 IRS 
5) that has a bright [Ne II] line detected in the Spitzer spec- 
trum, and an intermediate-mass Herbig Be star ( V892 Tau) be- 
cause its low resolution spectrum hinte d toward neon emission, 
but ju st below the detection threshold (Bal dovin-Saavedra et al.l 
1201 ll). The stars known to b e jet-drivi ng sources are: Co Ku Tau 
1 dEisloffel & Mundtlll998l). X Z Tau dKrist et alj|2008l). L155 1 
IRS 5 dRodrfguez et al.ll200l . and UY Aur dHktheTatl[T997l) . 
To our knowledge, the rest of the stars do not have reported out- 
flows or jets in the literature. In addition, among the targets se- 
lected the following are binari es not resolved by Spitzer. M HO- 1 
and MHO- 2 (3"9 separation. [Kraus & Hillenbrandll2009l) . V892 
Tau (0 705. ISrrdth et al.l 120051: iMonnier et al.ll2008l) FS Tau A 



1, jHa 

(0788, iMcCabeetal J 120061). XZ Tau (0730. lHaas et al] ["l990). 
V853 Oph (0'/3.lMcCabeetal] |2006l). and CoKu Tau 1 (0724, 
iPadgett et al.l fl999). TableQ]summarizes the stellar properties of 
the stars selected for this high-resolution spectroscopy follow- 
up. 

The scope of this study is to o btain high-resolutio n spectra 
of the [Ne II] line (12.81355 fim, lYamada et1d1[l985l) . and by 
determining the line center and studying the profiles, to obtain 
observational constraints on the emission mechanism of [Ne II]. 
The article is organized as follows: in Sect. 12. II we describe the 
VISIR observation strategy and data reduction, in Sect. 12.21 we 
present complementary observations obtained in the optical with 
UVES-VLT and its data reduction, in Sect. [3] we present the re- 
sults, the discussion in Sect. [4] and finally the conclusions in 
Sect. [5] 

2. Observations & Data Reduction 

2.1. VISIR-VLT 

Mid-infrared high-resolution s pectra were obtained with the im- 
ager and spectrometer VISIR (Lagage et al. 2004), installed at 
the Melipal telescope (UT3) of the VLT in three runs: May 
2009, January 2010, and January 2011. Spectra were obtained 
in high-resolution mode (HR) using the long-slit (3275 long). 
This configuration covers a spectral region between 12.793 and 
12.829 jum. One star (IRS 60) was observed using the cross- 
dispersed mode (HRX), that uses a short slit of 471, with the 
same spectral coverage. The slit width used was 074, achieving 
a resolving power of ~ 30000 in both configurations, measured 
from the FWHM of sky lines present in the spectra, which in 
velocity resolution corresponds to ~ 10 km s . The slit was 
oriented in the default North-South orientation, except for the 
sources known to be binaries with a separation such that the sin- 
gle components could be separated with VISIR. In those cases 
the slit orientation was adapted in order to obtain the spectra of 
the two components of the system. Standard chop-nodding along 
the slit was applied to correct for mid-infrared background emis- 
sion with a cho p throw between 8" and 12". Standard giant stars 
(from the list of ICohen et al.l[l999b or an asteroid (Psyche) were 
observed immediately before or after the science target to cor- 
rect for telluric absorption and to obtain the flux calibration^] A 



(0724, Hartigan & Kenyon 2003; Hio ki et alj 1201 ll). UY Au r 



1 In principle, giant stars are not ideal telluric standard stars for 
medium and high-resolution infrared spectroscopy due to the presence 
of photospheric absorption lines, and asteroids are better telluric stan- 
dards. In our case, however, giant stars were sufficient since Ne II detec- 
tions (before telluric correction) were very strong, as shown in Fig. IA.ll 
Photospheric absorption in giant star spectra are faint, after checking 
with the asteroid spectra. The detections are, therefore, not due to inap- 
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Table 1. Stellar properties of the sample studied 



Name 


a (hms) 


g (= < ») 


Distance " 


Spectral Type 


Ref. 


T 6 




(J2000) 


(J2000) 


(pc) 






(erg s _1 ) 


MHO 1 
MHO 2 


04:14:26.0 
04:14:26.4 


+28:06:04.2 
+28:05:59.7 


140 
140 


M2.5 
M2.5 


1 
1 


1.6 x 10 30 


V892 Tau 


04:18:40.7 


+28:19:12.3 


140 


B9 


1 


9.2 x 10 30 


CoKu Tau 1 


04:18:51.6 


+28:20:26.5 


140 


M0 


1 




FS Tau A 


04:22:02.2 


+26:57:30.2 


140 


M0 


2 


3.2 x 10 30 


SST042936+243555 


04:29:36.2 


+24:35:52.5 


140 


M3 


2 




AZ lau 


1)4:31: iy.a 


i 10.11. C7 -y 

+ 1o:1j:5/.j 


i yin 
14U 


\ /to 

M2 


3 


y.o x it) 


L1551 IRS 5 


04:31:34.2 


+ 18:08:04.7 


140 


<K6 


1 


2.1 x 10 28 


UY Aur 


04:51:47.9 


+30:47:23.0 


140 


M0 


3 


2.6 x 10 29 


VZCha 


11:09:24.4 


-76:23:22.4 


178 


K7 


4 


2.5 x 10 29 


RXJ1 11 1.7-7620 


11:11:47.0 


-76:20:10.8 


163 


K6 


5 


1.5 x 10 30 


V853 Oph 


16:28:45.3 


-24:28:18.9 


125 


M3.75 


6 


3.1 x 10 30 


IRS 60 


16:31:31.0 


-24:24:40.0 


125 


K2 


4 


2.6 x 10 29 


EC 92 


18:29:57.8 


+01:12:53.1 


260 


K7-M2 


7 


9.5 x 10 30 



" Distance com piled inlGudel et al.1 fcOlOh 

* L x come from Gii del et al.l J2007T) or lGudel et all d2010h 
Ref erences: mlLuhman et alj I201CI); (2)lRebull et al] d2010h : (3) lHartigan & Kenvonl i2003h : (4l lTorres etafl d2006h : (5) lLuhmarJ d2004r> : (6) 
IWilking et al] d2005l) ; (71 lPreibiscbl ( H"999h 



summary of the VISIR observations is presented in Table [2] in- 
cluding the date of the observations, the exposure time, the mode 
(HR high-resolution long-slit or HRX high-resolution cross- 
dispersed), the airmass at the beginning and end of the observa- 
tions, the seeing in the mid-infrared measured from the FWHM 
of the continuum, the position angle (following the standard con- 
vention of positive angles in the North-East direction), the star 
or asteroid used as calibrator, its respective exposure time, and 
airmass. 

The VISIR raw frames were processed with the standard 
VISIR data reduction pipeline version 3.7.2, using the command 
line application esorex. The output of the pipeline is a series of 
files including the extracted spectrum, a 2D image spectrum, 
a pixel-to-wavelength map, and a synthetic model spectrum of 
the calibration star. Starting f rom the 2D spectrum we have per- 
formed further steps based on lBoersma et al.l(l2009l) for the spec- 
tra obtained in HR mode. 

i) Any remaining background emission is corrected by taking 
each pixel row and fitting a second order polynomial, ignoring 
the source power spread function (PSF) profile. The result from 
the fit is then subtracted from each pixel row, obtaining a 2D 
spectrum that is background-corrected. 

ii) A weight map is created by collapsing the image in the 
dispersion direction, normalizing it, and then expanding it in the 
spatial direction. The science frame is then multiplied by the 
weight map. The spectrum is obtained by collapsing the result- 
ing image in the spatial direction. 

Hi) The wavelength calibration is performed using the atmo- 
spheric emission lines. This is done by taking a raw frame and 
cross-correlating it with a model of the atmospheric spectrumQ. 
Gaussian profiles are fitted to the atmospheric lines in order to 
obtain accurately the position of their center. By fitting a sec- 
ond order polynomial the pixel-to-wavelength conversion is ob- 
tained. This method gives a precision of 0.4 - 0.6 km s in the 
wavelength solution. 



propriate telluric correction. Nevertheless, we emphasize that observa- 
tions to detect faint Ne II features should ideally be done with asteroids. 

2 Available on the VISIR Tools website, see 
http : //www . eso . org/sci/facilities/paranal/instruments 



iv) The observed calibrator spectrum is corrected for differ- 
ences in airmass and air pressu re with the science targ et, follow- 
ing the procedure described in Carmona et al] (1201 lb . The cali- 
brator synthetic spectrum is divided by the observed calibrator 
spectrum. This result is then multiplied to the extracted spec- 
trum of the science target to correct for telluric absorption and 
to get an absolute flux calibration. The uncertainty in the flux 
calibration is of the order of 20 - 30%. 

Steps i, ii, and Hi were applied to both the science target and 
its calibrator. The VISIR spectra were corrected for barycentric 
and radial velocity. 

Whenever an asteroid was observed to correct for tel- 
luric absorption, the absolute flux calibration was based in 
Spitzer fluxes. In the particular case of UY Aur, the Spitzer flux 
corresponds to the binary, therefore to obtain an absolute flux 
calibration for the single components we used in addition the 
flux ratio between the c omponents, derived f rom high-resolution 
imaging in the N band dMcCabe et al]|2006l) . 

For the binaries, we converted the separation of the spectra 
from pixels to distance in the sky using the VISIR pixel scale 
(07127/pix). The binary systems resolved by our VISIR obser- 
vations are MHO-1/2, and UY Aur, for which the projected sep- 
arations observed are consistent with the separations reported 
in the literature. The small separation of XZ Tau did not allow 
the system to be spatially resolved, although the observed cross- 
profile was broad. The CoKu Tau 1 and FS Tau A binaries are 
also not resolved and show cross dispersion profiles consistent 
with a single source. Observations of EC 92 include the nearby 
star EC 95, located at a distance of 5", however the spectrum 
of EC 95 has a very low signal-to-noise ratio and it is not ex- 
ploitable. 

2.2. UVES-VLT 

After the end of our last VISIR run we requested observing time 
to obtain optical spectra with UVES-VLT for a number of targets 
with detected [Ne II] emission in VISIR, under DDT program 
286.C-5038(A). The purpose was to derive radial velocities at 
high precision, only possible from high-resolution spectra (see 
details in IA.2.U . Radial velocities are obtained by measuring 
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Table 2. Summary of the VISIR observations 



Name 


Date 


ton- source 


Mode 


Airmass 


MIR seeing 


Pos. Angle 


Calibrator (Sp. T) 


teal 


Airmass 






(s) 




start/end 


(") 


o 




(s) 




MHO 1/2 


2011-01-21 


1800 


HR 


1.71/1.89 


0.39 


-153.8 


HD 23319 (K2.5III) 


240 


1.02 




2011-01-23 


1800 


HR 


1.68/1.81 


0.33 


-153.8 


HD 29291 (G8III) 


240 


1.01 


V892 Tau 


2010-01-06 


3600 


HR 


1.71/1.70 


0.36 





HD 20644 (K4III) 


240 


1.69 


CoKu Tau 1 


2010-01-04 


3600 


HR 


1.72/1.69 


0.32 





HD 17361 (K1.5III) 


240 


1.70 


FS Tau A 


2010-01-04 


2760 


HR 


1.73/2.14 


0.33 





HD 27639 (M0III) 


240 


1.95 




2010-01-05 


3600 


HR 


1.71/1.62 


0.29 





HD 27482 (K5III) 


240 


1.70 


SST042936+243555 


2011-01-21 


1800 


HR 


1.78/2.17 


0.39 





HD 27639 (M0III) 


600 


2.46 


XZTau 


2011-01-22 


1800 


HR 


1.52/1.78 


0.32 


-138.6 


Psyche 


120 


1.85 


L1551 IRS 5 


2010-01-05 


3600 


HR 


1.41/1.78 


0.28 





HD 27482 (K5III) 


240 


1.64 


UY Aur 


2011-01-22 


1800 


HR 


1.76/1.81 


0.36 


132.4 


HD 29291 (G8III) 


240 


1.00 




2011-01-23 


1800 


HR 


1.88/2.14 


0.38 


132.4 


Psyche 


120 


1.77 


VZ Cha 


2011-01-22 


1800 


HR 


1.87/1.76 


0.32 





Psyche 


120 


1.85 


RXJ1 11 1.7-7620 


2011-01-21 


960 


HR 


1.79/1.74 


0.39 





HD 27639 (M0III) 


600 


2.46 


V853 Oph 


2009-06-01 


6120 


HR 


1.01/1.50 


0.27 





HIP 90185 (B9.5III) 


600 


1.11 


IRS 60 


2009-05-31 


3600 


HRX 


1.19/1.01 


0.40 





HIP 84012 (A2IV) 


840 


1.05 




2009-06-01 


1800 


HRX 


1.06/1.00 


0.40 





HIP 84012 (A2IV) 


480 


1.04 


EC 92 


2009-05-31 


4080 


HR 


1.15/1.18 


0.26 


-173.7 


HIP 99473 (B9.5III) 


1200 


1.15 




2009-06-01 


3600 


HR 


1.36/2.24 


0.32 


-173.7 


HIP 84012 (A2IV) 


840 


1.38 



Table 3. Summary of the UVES observations 



Name 


Date 


Exp. time 


Airmass 


Seeing 






(s) 




(") 


CoKu Tau 1 


2011-03-02 


900 


2.0 


1.4/1.2 


FS Tau A 


2011-03-02 


900 


2.1 


1.2/1.3 


V853 Oph 


2011-03-02 


300 


1.3 


1.3/1.1 



the shift of the absorption lines of the stellar atmosphere with 
respect to a reference spectrum and are crucial to determine the 
rest velocities of the [Ne II] lines detected with VISIR. UVES 
(Ultraviolet and Visual Echelle Spectrograph) is located at the 
Nasmyth platform B of the second Unit Telescope (Kueyen) of 
the VLT, offers cross-dispersed echelle spectra between 300 up 
to 1100 nm. The light beam is splitted into two arms; UV-Blue 
and Visual-Red. These two arms can operate separately or in 
parallel. We used the red arm covering the wavelength range 
500 - 595 and 605 - 700 nm. This setting allowed us to cover 
a region rich in absorption lines, required for the determination 
of the stellar radial velocity, and including important emission 
lines, such as [O I] and Ha. The slit width used was 1'.'2, achiev- 
ing a resolving power R ~ 33000. The summary of the observa- 
tions is presented in Table [3] 

We reduced the data following the standard UVES pipeline 
recipes, version 4.9.0, using the command line application es- 
orex. Optical data need to be corrected for several instrumental 
effects. To correct for the electronic noise of the camera and pos- 
sible systematics, a short exposure (bias) is taken with the shut- 
ter closed; the bias frames give the read out of the CCD detector 
for zero integration time. Several bias frames are combined into 
a master bias; in the case of UVES five bias frames are taken. 
The more bias exposures are taken, the less noise will be intro- 
duced into the corrected images. Furthermore, as the telescope 
does not illuminate the detector homogeneously, and the quan- 
tum efficiency of the CCD is not necessarily the same over all the 
pixels, an exposure of an homogeneously illuminated area (flat 
field) needs to be taken. Five flat fields are then combined to cre- 
ate the normalized master flat field for UVES observations. The 
bias frame is subtracted of each raw science frame and the resul- 
tant spectrum is divided by the master flat field frame. In gen- 



eral, a correction for the dark current of the CCD is also needed, 
and this is done by taking a long exposure with the shutter closed 
(dark frame). In the case of UVES, the detector dark currents can 
be considered negligible and be excluded from the data reduction 
process. The wavelength calibration of UVES spectra is done by 
using observations of a ThAr lamp and is part of the standard 
procedures of the UVES pipeline. The error in the determina- 
tion of the wavelength solution is 2 mA, typical of the UVES 
pipeline. The spectra were also corrected for radial and barycen- 
tric velocity. Given that we did not observe standard stars, the 
correction for the telluric absorption features and absolute flux 
calibration was not possible. However, the data reduction pro- 
cess includes the correction of the spectra for the sky airglow 
emission. We used the R magnitude and Johnson zero point pho- 
tometry to determine the flux density in the continuum. 

3. Results 

3.1. [Ne II] detections in VISIR spectra 

The [Ne II] emission line was detected in seven stars in 
our sample: MHO-1, V892 Tau, CoKu Tau 1, FS Tau A, 
SST 042936+243555, V853 Oph, and EC 92. The line was fitted 
using a Gaussian profile plus a linear component to account for 
the continuum. All the parameters of the function were left free 
to vary. For the cases with non detections, an estimation of the 
upper limit was calculated as three times the standard deviation 
of the continuum flux (cr) times the spectrograph full width half 
maximum (FWHM). Figure [TJ displays the spectra after wave- 
length and flux calibration, plotted including the error bars in 
the flux. The green line represents the Gaussian fit, and the red 
dot-dashed line is the 3cr upper limit threshold. The spectra are 
plotted in the stellocentric frame; the vertical dashed line shows 
the position the lines would have if centered at the stellar ra- 
dial velocity. The figure was divided in two panels: detections 
in the upper panel and non detections in the lower one. In addi- 
tion, the detections are separated in two groups: stars showing 
a line blueshifted to high velocities, and stars with a line nearly 
centered or shifted by a few km s . Figure lATI displays the 
observed spectra of the stars and its telluric standard before ap- 
plying any velocity correction. Only one data set is included for 
the stars observed in two different days. 
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Fig. 1. Spectra of the stars observed with VISIR plotted in stellocentric frame and rebinned by a factor of two (average spectra are 
shown when there are more than one observations). The detections are plotted in the upper panel, non detections in the lower one. 
The green line is the Gaussian fit to the line, the red dot-dashed line is the 3cr detection threshold per pixel. The error bars shown 
are 1<t values per pixel. 



Table [4] shows the results obtained from the fits to the 
[Ne II] emission line: the integrated line fluxes or upper limits 
in Col. 2, the FWHM of the detected lines in Col. 4, and the cen- 
ter of the Gaussian in Col. 5. In addition, we report the Spitzer 
fluxes for comparison in Col. 3, the stellar radial velocity and its 
reference in Cols. 6 and 7, respectively. The errors of the mea- 
surements are reported between parentheses. 

The detection of [Ne II] from the Herbig Be star V892 
Tau that we present here is the first report of [Ne II] emis- 



sion from an intermediate mass pre-main sequence star. The 
[Ne II] line is detected close to the stellar rest velocity in three 
cases: CoKu Tau 1, FS Tau A, and V892 Tau. For them, we stud- 
ied the sy mmetry of the line pro files by using the procedure de- 
scribed in Pasc ucci et al.1 d201 ll) : each line is flipped and shifted 
in wavelength direction to get a good match with the original 
profile. The two profiles are then subtracted and the result is 
plotted; the more symmetric the line, the closer to zero the dif- 
ference. The results are shown in Figure |2] in the upper panel 
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Table 4. Results for the [Ne II] line obtained from the VISIR spectra. The peak is given in the stellocentric frame. Errors are reported 
between parentheses. 



[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


Name 


Flux 


Spitzer Flux * 


FWHM 


Peak 1 ' 


vrad 


vrad Ref. 




(10 -15 ergcnT 2 s _1 ) 


(10~ 15 erg cirT 2 s~ 


- 1 ) (kms- 1 ) 


(km s~') 


(km s _1 ) 




Detected lines 


MHO 1 


8.0(1.5) 


120* 


26.7 (3.3) 


-121.8 (6.6) 


16.0 (6.4)* 


Bertout & Genova (2006) 


V892 Tau 


179 (16) 


< 180 


26.0 (3.4) 


-2.13(6.5) 


16.0 (6.4)* 


Bertout & Genova (2006) 


CoKu Tau 1 


20.4 (2.5) 


120 


55.2(3.3) 


3.6(1.5) 


15.0 (0.8) 


White & Hillenbrand (2004) 


FS Tau A 


16.6(1.7) 


77 


26.8 (1.7) 


-2.9 (0.7) 


17.1 (0.1) 


This work 


SST042936+243555 


7.8 (1.7) 


11 


26.0 (3.8) 


-89.1 (6.6) 


16.0 (6.4)* 


Bertout & Genova (2006) 


V853 Oph 


3.7 (0.8) 


16 


26.5 (3.9) 


-35.8 (1.6) 


-8.9 (0.1) 


This work 


EC 92 


4.2(1.2) 


15 


16.2 (4.3) 


-31.3 (2.3) 


-7.1 (1.5) 


Covev et al. (2006) 


Upper limits 


MHO 2 


< 0.5 


120* 






16.0 (6.4)* 


Bertout & Genova (2006) 


XZTau 


< 5.2 


32 






17(7) 


Folha & Emerson (2000) 


L1551 IRS 5 


< 4.3 


580 






8.2(1.5) 


Covev et al. (2006) 


UY Aur A 


< 2.4 


49* 






18(3) 


Barbier-Brossat & Figon (1999) 


UYAurB 


< 1.2 


49* 






18(3) 


Barbier-Brossat & Fison (1999) 


VZ Cha 


< 2.3 


3.8 






19.0(1.7) 


Torres et al. (2006) 


RXJ1 11 1.7-7620 


< 3.3 


5.1 






15.5 (0.3) 


Torres et al. (2006) 


IRS 60 


< 1.6 


14 






-2.2 


Torres et al. (2006) 



* Spitzer fluxes come from lGudel eTai] d2010h or lBaldovin-Saavedra et all fcOllh 

t The error in the center was calculated considering the contribution of the error in the Gaussian fit and in the stellar radial velocity 

* In these cases there is no radial velocity in the literature, the average radial velocity of Taurus was assumed. 

The binary systems UY Aur (A,B) and MHO 1/2 cannot be separated in the Spitzer spectrum, therefore the flux reported is for the system. 



of the figure, the original and the flipped spectrum are shown, 
while the lower panel displays the difference between both pro- 
files. This method does not show any evident asymmetry in the 
line profiles. 



The spatial extent of the emission can be quantified by sub- 
tracting the average PSF observed in the continuum from the 
2D spectrum, as presented in the position velocity diagrams in 
Figure [3] (contours start at 3 <x with an increase of 5<x at each 
step), bearing in mind the limited spatial resolution imposed by 
the instrumental PSF. CoKu Tau 1 presents symmetric emission 
centered at the stellar rest velocity, with no extended emission 
beyond the (X'3 (40 AU at the distance of Taurus) spatial reso- 
lution of the VLT in the mid-infrared. The average PSF in the 
continuum is also co nsistent with a poin t source, indicating that 
the (X'24 companion dPadgett et al.l fl999) was not detected. The 
centroid of the emission at negative velocities is slightly dis- 
placed to "negative" spatial extension (about 10 AU) whereas 
the centroid at positive velocities goes to "positive" spatial ex- 
tension (also about 10 AU). This could be a signature of disk 
rotation, but considering the fact that the slit was oriented with 
PA=0°, i.e., close to the direction of the bipolar jet (PA=210°), 
the centroid's displacement could also reflect emission by the 
bipolar jet. We address this point in the discussion (Section l4~T1 ). 
The emission from FS Tau A is centered at negative velocities, 
with evidence of extension toward negative velocities, indica- 
tion of photoevaporative wind. The symmetry of the emission 
is much less evident than in CoKu Tau 1 . The 3cr level extends 
up to velocities of +40 km s . The [Ne II] PV diagram in FS 
Tau A suggests some extended emission be yond Q'/3 that cannot 
be att ributed to the close (J.'24 companion (Hartigan & Kenvorl 
120031) . The [Ne II] emission from V892 Tau shows an extension 
toward the blue side of the spectrum. Its PV diagram also indi- 
cates extended emission beyond Of.' 3. 



3.1 .1 . Comparison with Spitzer line fluxes 

We now compare the line fluxes obtained with VISIR with the 
fluxes obtained previously from Spitzer low-resolution spectra. 
In general, VISIR fluxes tend to be lower than Spitzer fluxes 
(reported in Table |4] column 3). Focusing only on the detected 
lines, we found that in one case (V892 Tau) the flux obtained 
with VISIR corresponds to the Spitzer upper limit. The differ- 
ence between Spitzer and VISIR fluxes can go as high as one 
order of magnitude. In the case of MHO 1/2, the Spitzer flux 
reported is for the unresolved binary, but the [Ne II] line was de- 
tected only in MHO 1 with VISIR. We found that for FS Tau A 
and SST042936+243555 the VISIR flux is between 30 and 40% 
lower than the Spitzer flux. For CoKu Tau 1, V853 Oph, and 
EC 92 the difference between Spitzer and VISIR fluxes are high, 
the flux obtained with VISIR observations represent between 
20 - 30% of the Spitzer flux. The variation b etween Spitzer 
and V ISIR fluxes has been previously reported bv lPascucci et alJ 
(1201 lh. where the y found on the order of 30% for TW Hya. 
Sacc o et alJ d2012l) also found evidence of the flux discrepan- 
cies in their targets for Class I and Class II stars, arguing that 
the [Ne II] emission is extended at least in Class I stars, and 
partially in Class II stars. However, they found flux ratios be- 
tween Spitzer and VISIR that are consistent within a factor of 2 
for pre-transition and transition disks, suggesting that the emis- 
sion is produced close to the star, within ~ 20 - 40 AU. In our 
study, the reason for the observed variation remains unclear, al- 
though plausible but unlikely explanations might be unidenti- 
fied flux calibration issues, slit losses in VISIR observations due 
to an incorrect centering of a star in the slit, unexpected tele- 
scope drifts, or to atmospheric conditions. A faint, very broad 
line could also remain undetected. The most likely explanation 
for the flux discrepancies is the different beam sizes between 
Spitzer and VISIR that translate into possible extended emission 
unresolved and detected in Spitzer but undetected with VISIR. 
Indeed, Spitzer spectra are extracted from a region covering the 
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Fig. 2. Upper panel: Plotted as grey filled area the profile of the [Ne II] line and in red solid line the profile flipped. Lower panel: 
Difference between the original and the flipped profile. 
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Fig. 3. Upper panels: Observed VISIR spectrum of CoKu Tau 1 , FS Tau A, and V892 Tau. Lower panels: Position velocity diagrams 
after subtraction of the continuum PSF. The contour levels start at 3<x and increase with a step of 5cr. The PSF FWHM is * (X'3. 



whole slit width, for IRS high-resolution this means a region of 
-11", i.e., about 1500 AU at 140 pc. 

We note that L1551 IRS5 was selected for this program be- 
cause a bright [Ne II] line was detected in the Spitzer spectrum, 
with a flux of 5.8xl0~ 13 erg s _1 crrT 2 , i.e., a luminosity of the or- 
der of L[Ne ii] ~ 10 3l) erg s _I . Although the line was not detected 
in spectroscopy mode, the target was detected in the [Ne II] nar- 
row band image, likely due to the strong continuum emission. 
This Class I star is known to drive a highly collimated bipo- 
lar outflow observed in the optical, near- infrared and radio (e.g., 
iDavis et alJl200llPvo et ai]|2009l:IWu etai]|2009l) . The position 
angle of the outflow is measured at 260° fo r the northern co mpo- 
nent and 235° for the southern component dPvo et al.l l2009). The 
VISIR observations were performed with a P.A.= 0°. The non- 
detection of [Ne II] in the spectrum suggests that the emission 



originates in the outflow, at distances larger than « ff.'4 from the 
central star, although the short integration in the acquisition im- 
age did not allow the detection of extended emission. We did not 
repeat the observations at a different PA. because the outflow 
explanation is the most likely in this case. 

3.2. Optical emission lines detected in UVES spectra 

Low-excitation forbidden lines are used as diagnostics 
of the physical conditions of the gas in the wind and 



shock environments of T Tauri s tars (e.g. Hamariiil 1994; 
iGomez de Castro & Pudritzl Il993l lOuve d & Pudritzl 11994 . 
For example, the ratio [S II]^6716/[S II]^6730 is strongly 
dependent on the electron density n e , while the ratios 
[SII]/16730/[OI]/16300 and [N II]/16584/[0 1]/16300 are 
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dependent on the electron temperature, T e . The [O II /16300 lin e 
is usually observed in T Tauri spectra (e.g., ICabrit et"aT1 ll990). 
often exhibiting a broad high velocity component blueshifted up 
to several hundred km s _1 (HVC) and a low-velocity component 
(LVC) w hich is narrower and blueshifted by only a few km s _I 
(see e.g jHartigan et al.lll995l) . The LVC and HVC are thought 
to have different origins; the HVC is believed to be tracing a jet, 
while the LVC is interpreted as a tracer of a photoevaporative 
disk wind by lErcolano & Owenl (l2010h . In particular, these 
authors proposed the presence of both [Ne II] and [O I] as an 
indicator of an X-ray photoevaporative disk wind. Permitted 
emission lines of Fe II are als o observed in the optical. For 
example, Beri stain et al.l (Il998l) found line profiles showing a 
narrow and a broad component that were interpreted as being 
produced by different kinematic zones: the narrow component 
is produced by turbulence and broadened by stellar rotation, 
originating in gas located in the post-shock region, and the 
broad component is produced in gas infalling in the accretion 
funnel. Furthermore, He I emission lines are presented as a 
powerful diagnostics, because their high excitation potential the 
emission is restricte d to regions of high t emperature or close to 
the ionizing source (Beristain et al. 20QTJ). 

The optical spectra obtained with UVES show several of 
these emission lines, whose profiles can be compared with the 
profile of the [Ne II] line from VISIR to understand the circum- 
stellar environment of the sources and provide hints on the origin 
of the [Ne II] lines observed. In Table [5] we present the proper- 
ties of the optical lines detected in the spectra of CoKu Tau 1, 
FS Tau A, and V853 Oph. In the second column we include the 
rest wavelength of the lines taken from the NIST Atomic Spectra 
Database Q. The line center with respect to the systemic velocity 
of the star, and FWHM (Cols. 3 and 4) were obtained by fit- 
ting a Gaussian profile to the lines, and the errors reported come 
from the fit. The line fluxes (Col. 5) were calculated using as 
continuum the stellar flux derived from the R magnitude and the 
Johnson zero point photometry. In Col. 6 we present the equiva- 
lent widths (EW) calculated by integrating the line fluxes within 
3<r of the Gaussian fit. 



4. Discussion 

4.1. Origin of[Ne II] emission: observations and models 

The emission of [Ne II] from protoplanetary disks has been 
discussed in a series of articles that model the irradiation of 
the disk by high energy photons, each of them gives a differ- 
ent possible scenario for the formation of this emission line 
(seeblassgold et al.l2007tlAlexanderl2008tlGorti & Hollenbachl 
2008; Ercolan o & Owenl |2Q1Q|> . Observations at high spectral 
and spatial resolution provide line profiles and shifts that can 
be used to test the d ifferent models. 

Glassgold et al.l d2007l) predicted [Ne II] emission from an 
X-ray irradiated disk atmosphere. In this model, the [Ne II] emis- 
sion comes from a region within 20 AU from the star. The pre- 
dicted line profile is double peaked with extended wings asso- 
ciated with material at small radii, although the authors did not 
giv e any estim a tes of the line widths expected. 

[Alexander (2008) modeled the [Ne II] line profiles from 
a photoevaporative disk wind driven by stellar EUV photons, 
assuming a spectral resolution of R ~ 30000. The model in- 
cluded different disk inclination angles. The line profiles ob- 
tained varied with the disk inclination considered. A disk viewed 



3 http://physics.nist.gov/ 



face-on (inclination angle i = 0°) showed a single peaked line 
blueshifted by typically < 10 km s . On the other hand, a disk 
viewed edge-on (i = 90°) showed a double-peaked profile. The 
emission is expected to be centered at the stellar velocity with 
one peak blueshifted and the other redshifted. The line widths 
predicted are FWHM ~ 10 km s -1 for the face-on case and 
FWHM ~ 30 km s for the edge-on case. In addition, a disk 
with an inner hole of 9 AU (in both dust and gas) is also mod- 
eled. For the face-on case, the line width and shift are equiva- 
lent to the disk without an inner hole. In the edge-on case, the 
line was narrower for the model with an inner hole than for the 
model without an inner hole (FWHM ~ 20 km s -1 ). In addi- 
tion, the line was double peaked and centered at the stellar rest 
velocity. A narrower line is in fact expected since the material 
closer to the star has larger rotational velocity and gives rise to 
the broadening of the l i ne. 

Ercolano & O wenl (1201 Ol) considered the photoevaporation 
of the disk driven by EUV and X-rays and modeled the pro- 
files for a series of emission lines, putting special attention on 
[Ne II] and [O I] (/16300). The model considered a disk with and 
without inner holes of 9, 14, and 30 AU at different inclination 
angles between 0° and 90°. Furthermore, the X-ray luminosity 
was also variable between log(Lx) = 28.3 to log(Lx) = 30.3. 
The profile of the [Ne II] line for the model without an inner hole 
varied according to the disk inclination. A face-on disk showed 
a narrow profile, with a FWHM between ~ 4 km s _1 for low 
Lx and ~ 10 km s _1 for high Lx- For the edge-on case, the line 
profiles are larger but do not change much with Lx', the FWHM 
are between ~ 20 km s _1 and ~ 23 km s _1 for the low and 
high Lx, respectively. Only in the case of a face-on disk and 
low Lx was the line found to be centered at the stellar velocity; 
in all other cases the line was blueshifted by up to ~ 5 km s _1 . 
All the profiles obtained were asymmetric, with a shoulder to- 
ward the blueshifted part, which became more evident for the 
high Lx model. 

Using the predictions of such models we can attempt to give 
an interpretation to the line profiles observed. Following the or- 
der of the panels of Fig. [T]we can identify three regimes for the 
[Ne II] emission: 

Highly blueshifted [Ne II] emission from jets: The line pro- 
files shown in the upper panel of Fig. Q]correspond to this group 
and they are characterized by a large velocity shift toward the 
blue. The emission is likely to originate mainly from shocks in 
jets. For MHO-1 and SST 042936+243555 there is no available 
measurement of the radial velocity (we assumed the average ra- 
dial velocity of the Taurus Molecular Cloud) but since they are 
observed with high velocity shifts, -122 and -89 km s _1 respec- 
tively, we conclude that the emission is likely originating in a jet. 
V853 Oph and EC 92 also likely belong to this group, with ve- 
locity shifts of -36 and -31 km s _1 respectively, despite the low 
peak velocities, perhaps due to inclination effects. 

There is no record in the literature of jets detected in our 
sample of stars with highly blueshifted [Ne II] emission. In the 
case of SST 042936+243555, this is likely due to the low num- 
ber of observations dedicated to look for jet or outflow emission 
published. From a search in the literature, spectroscopic obser- 
vations in the optical or in the near-infrared do not show any 
evidence of outflow or jet emission in MHO-1, V853 Oph, or 
EC 92. 

IShang et ail d2010l) modeled the profiles of [Ne II] and 
[O I] (/16300) expected from a jet. The centroid of the lines show 
a large shift toward the blue, reaching -200 km s , for an in- 
clination angle of 45°. The lines are asymmetric with a broad 
shoulder toward the red side of the spectrum. The excess emis- 
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sion toward the red can reach up to +100 km s _1 , depending on 
the inclination angle. The profiles of the model were calculated 
for a resolution of 1 km s _1 , much higher than the resolution 
achieved with our VISIR observations. However, the dominant 
feature is a large blue-shifted peak, like observed in MHO-1, 
SST 042936+243555, V853 Oph, and EC 92. 

Weakly blueshifted [Ne II] emission from photoevaporative 
wind: The centroids of the lines in this category are close to the 
stellar rest velocity or show a small velocity shift. To this group 
belong FS Tau A and V892 Tau. In particular, the [Ne II] pro- 
file from V892 Tau shows an excess toward the blue, and a line 
width of 25 km s _1 . The star is known to be a binar y of separation 
~ 5 A U. The circumbinary disk was detected by Mo nnier et al.l 
(2008) through infrared imaging. The inclination of the disk is 
estimated to be i = 60° and the disk inner hole to be ~ 35 AU. 
FS Tau A is also a close binary system, with separation 0V24 
dHartigan & Kenvon|[2003l equivalent to 34 AU at the distance 
of Taurus) and surrounded by a c ircumbinary disk that extends 
up to 630 AU dHioki et alj|201 ll) . The spatial extension of the 
emission (Figure in both systems shows an excess toward the 
blue side of the spectrum, i.e., an indication of a photoevapora- 
tive disk wind. 

The special case of CoKu Tau 1: The emission line de- 
tected in CoKu Tau 1 shows a flat-topped, broad profile with 
a FWHM of 55 km s -1 , and wings extending up to velocities 
close to ±80 km s _1 . The emission is also spatially unresolved. 
The immediate interpretation is that the [Ne II] line originates 
from the disk atmosphere, from material located at small disk 
radii. H owever, the star is kno wn to be viewed close to edge-on 
(i = 87°, Robit aille et alj2007l) and is known to have a bipolar jet 
detec ted at both positive and negative velocities relative to th e 
star dMovsesvan & Magakvan|[l989t lEisloffel & Mundd 1 1998h . 
Photoevaporation models at nearly edge-on inclinations predict 
large wings, but due to mo del construction they never extend 
beyon d +40 km s -1 (e.g., lAlexanderl 120081: lErcolano & Owed 
120101) . Perhaps a fast photoevaporativ e wind could exp lain the 
broader wings observed in CoKu Tau 1 . Bas t et all (1201 ll) argued 
also that broad centrally peaked CO line profiles in protoplane- 
tary disks could not be explained by inclined disks in Keplerian 
rotation, but that a combination of emission from the inner part 
of the disk and a slow moving disk wind could explain the line 
profiles. Although the disk interpretation for CoKu Tau 1 is more 
physically interesting, the bipolar jet explanation is as valid: first, 
its [Ne II] flux is rather high with Spitzer and VISIR and is more 
typical of fluxes detected from jets. Second, although the PV 
diagram shows spatially unresolved emission, a careful analy- 
sis indicates peaks at "negative" spatial extension for negative 
velocities, and peaks at "positive" spatial extension for positive 
velocities. Considering the fact that we used a NS slit orientation 
with VISIR (and UVES) while the jet is oriented with PA=210° 
with positi ve velocities toward the South (and inversely for the 
counterjet; lEisloffel & Mundtll 19981) . the centroids of the [Ne II] 
peaks are consistent with the bipolar jet interpretation. On the 
other hand the direction of rotation of the disk in CoKu Tau 1 is 
unknown, and the [Ne II] peaks could also trace this rotation. In 
fact, the peaks extend up to ±10 AU. In conclusion, the origin of 
[Ne II] emission in CoKu Tau 1 is difficult to ascertain and could 
be either from the disk atmosphere or from the bipolar jet, with a 
possible contribution by a photoevaporative wind. A comparison 
with the optical forbidden line (see below) does not allow us to 
constrain better the origin of the [Ne II] line. 



4.2. Comparison between [Ne II] profile and optical lines 

We can further compare the profiles of the [Ne II] line with the 
profiles of the optical lines for the three targets observed in our 
UVES program. Figures |4] [5] and|6]show in black line and grey 
area the profile of the [Ne II] line and in red the profile of the 
optical lines detected. The spectra were continuum subtracted 
and normalized to the peak of the line to allow comparison. 

- CoKu Tau 1 

The profiles of the optical forbidden lines are generally com- 
patible with the profile of the [Ne II] line, i.e., similar center 
and FWHM. The [N II] lines are slightly broader than [Ne II], 
in contrast to the [O I] and [S II] lines which are of similar 
width (or slightly smaller in the case of oxygen). The Ha line 
is much broader, possibly tracing additional regions than the 
optical forbidden lines. The Fe II line (/16432) is also much 
broader (FWHM ~ 250 km s ). The optical forbidden lines 
are only shown in Figure |4] whereas the Fe II and Ha lines ap- 
pear in the Appendix (Fig. IA.2l i. The shapes of the two detected 
[N II] lines (/16548 and ,16583) are more difficult to interpret, 
they appear double or multiply peaked but this might also be 
due to a lower signal-to-noise ratio in that part of the spectrum. 
W e note that the optical l ine p rofiles are compatible with those 
of lWhite& Hillenbrand! J2004T) from HIRES data. The optical 
emission spectrum is completely dominated by optical forbid- 
den li nes (and Ha), which is typi cally seen in edge-on Class I 
stars dWhite & Hillenbrandl l2b04). The origin of these optical 
lines is unclear but could either be due to the disk atmosphere or 
from the inner bipolar jet. The profile of the [Ne II] line is con- 
sistent with the forbidden lines (Fig. [3}, which suggests that the 
line may either originate from the same region. In the bipolar jet 
interpretation, the broader [N II] lines may include a faster com- 
ponent, possibly because they would form closer to the faster jet 
axis, whereas the [O I], [S II], and [Ne II] lines form closer to 
the slower jet surface. The comparison of the optical forbidden 
lines and the [Ne II] line does not allow to determine conclu- 
sively the origin of the lines, although the bipolar jet origin is 
the most likely one. 

- FS Tau A 

The profile of the [Ne II] line in FS Tau A is generally similar 
to the profiles of the optical forbidden lines, and to He I (/15876; 
the line is, however, broader) and Fe II (^15 197; the line is red- 
shifted with respect to [Ne II] and the stellar radial velocity). The 
optical forbidden lines are nearly centered at the stellar radial 
velocity or slightly blueshifted with a FWHM ~ 30 km s _1 . The 
profiles are consiste nt with a photoevaporativ e wind flow, like 
the ones modeled by lErcolano & Owed (1201 Oh . nearly centered 
with a shoulder toward blueshifted velocities. This blueshifted 
shoulder is larger for [N II] (/16583). The observed line pro- 
files are consistent with the modeled profi les considering th e 
disk inclination, estimated at i = 30-40° dHioki et alj 1201 ll) . 
Furthermore, the presence of [Ne II] and [O I] are indica- 
tions of an X-ray driven photoevaporative flow, according to 
lErcolano & Owed (12010b . In fact, the low-velocity component 
of the [O I] fine can only be reproduced in the p resence of 
X-rays . A comparison with the profile modeled by lAlexanderl 
(2008) is more difficult since only the profiles for the edge-on 
and face-on disks are shown. Nevertheless the FWHM for the 
case i = 45° is reported and is compatible with the value de- 
rived from the observations: FWHM mo( j e i= 28 km s _1 versus 
FWHM obs = 27kms- 1 . 
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Fig. 4. Lines detected in the UVES spectrum of CoKu Tau 1 plotted as a filled grey area compared to the profile of [Ne II] (red line). 
The spectra were continuum subtracted and normalized to the peak of the lines to allow comparison. 



- V853 Oph 

V853 Oph shows a [Ne II] line blueshifted by 36 km s" 1 , 
with a FWHM — 26 km s -1 . The large velocity shift observed 
suggest that the line is emitted by a jet. In this star, the profiles 
of the optical lines are not always compatible with the profile of 
the [Ne II] line; we can classify them according to their shifts 
and widths: 

Lines emitted by the star or the disk: These are He I 
(,15015, 5876, 6678) and Fe II (/15018, 5169). They show differ- 
ent shapes; the He I lines are asymmetric with an excess toward 
redshifted emission, while the Fe II lines are well centered and 
narrow. 

Lines compatible with photoevaporation: These are [O I] 
(^5577,6300,6363). The [O I] lines show different shapes; the 
/15577 line is fairly well centered at the stellar rest velocity and 
shows the blue shoulder typical of photoevaporative wind flows, 
and is also compati ble with the disk inclination angle (z = 32°, 
lAndrews et al . 2010). The /16300 line, on the other hand, shows 
two velocity components: one close to the stellar rest velocity 
and the other with a higher velocity shift toward the blue, prob- 
ably due to a jet. The separation between the components are 
believed to be real, although there is likely a moderate absorp- 
tion between the components due to water molecular absorption. 
Indeed, although oxygen emission lines are present in the sky 
spectrum, they are removed in the data reduction process. But 
molecular absorption lines from the atmosphere can only be cor- 



rected with observations of a telluric standard. Since we do not 
have a telluric calibrator for the optical spectra, we used archival 
spectra of early type stars to look for contamination of atmo- 
spheric absorption lines in the profile of [O I]. There is, indeed, 
a 20-30 % absorption feature at the proper position, suggesting 
that some but not all of the absorption feature observed in the 
/16300 line is due to a water absorption line. 

Lines emitted in a jet: these are the [N II] (/16583) and the 
[S II] lines, and the high-velocity component of the [O I] /16300 
line. The [S II] (,16730) line matches the [Ne II] profile, with 
similar center and FWHM, whereas the [N II] line is slightly 
broader and shows emission at higher speeds, but the signal-to- 
noise ratio is not large. The l i ne pro files differ from the pro- 
files predicted by Shan g et aD (l2010h : they are narrow and the 
peak is located at velocities lower than predicted by the model; 
-37 km s" 1 against ~ 100 - 200 km s" 1 . 

In conclusion, at least in the case of V853 Oph, the detected 
[Ne II] line is not a tracer of photoevaporative wind since it traces 
a high-velocity component of a jet. Possibly the [Ne II] compo- 
nent from a photoevaporative wind exists, but it is much fainter 
than the jet component detected with Spitzer and VISIR. 

4.3. Summary of [Ne II] ground based observations. 



Putting together previous results (iHerczeg et ail 2007; 
I van Boekel et alj 120091 iNaiita et alj 120091; IPascucci & Sterzikl 
2009; IPascucci et alJ 1201 lb ISacco et aH 120121) with the ones 
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Fig. 5. Profile of the optical lines detected with UVES and the profile of the [Ne II] line for FS Tau A. The color code used is the 
same as Fig. [4] 



from our program, there is a sample of 26 objects, for which 
[Ne II] emission is detected from the ground with high spectral 
resolution to derive th e kinem atics of the [Ne II] emission. The 
study by ISacco et al] (|2012|) . published while this paper was 
being refereed, provided 12 new detections among 32 targets 
and was able to provide some conclusions about the origin 
of [Ne II] emission, in particular for Class I sources and for 
transition and pre-transition disk sources. Our study focusses 
more on Class II stars and provides additional comparison of 
the [Ne II] line emission with high-resolution optical lines, 
in particular forbidden lines such as [O I], [S II]. We can 
attempt to extract some general conclusions from the results 
we have in hand up to now, keeping in mind that ground based 
observations are challenging in terms of sensitivity. Studies 
based on Spitzer observations showed a weak correlation 
between the [Ne II] luminosity and the X-ray luminosity, in 
particular for optically thick disks sources. This correlation 



is expected if X-r ays are responsible for the neon emission. 
Sacc o et a D (120 12h found no obvious correlation between the 
[Ne II] and X-ray luminosities in their sample. Restricting to 
all detected sources in which ground based observations give a 
line consistent with disk emission and the ones consistent with 
photoevaporation, there is still no clear eviden ce of a trend. But, 
if dropping the 5 targets in lSacco et al. with an observed 

blueshift less than 18 km s _1 and whose disk inclinations are 
unknown, there is a possible trend for higher neon luminosities 
with higher X-ray luminosities. However, we emphasize that 
the sample is relatively small, and there are many uncertainties 
(classification, inclination, variability, etc). 

Another lesson from Spitzer observations is that bright neon 
emission is expected in jet sources, and that w henever jets 
are present, they dominate over the disk emission (Gude fet al.l 
201(4 iBaldovin-Saavedra et al.l l201lh . Ground-based observa- 
tions have shown that jets are significant emitters of [Ne II] line 
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same as Fig. [4] 



emission (e .g., [van BoekeJ_etal ] l2009l lP ascucci & Sterzik 2009, 
ISacco et al.ll2012i this work). The detected lines in jet sources 
tend to be fainter than lines compatible with d isk or slow wind 
emission (see also Fig. 4 in lSacco et al.ll2012h . This can be ex- 
plained because the larger beam of Spitzer was able to capture 
most of the extended emission produced by shocks, but this 
emission is filtered out by the narrow slit used in ground based 
observations. Jet activity is closely related to the evolution of 



pre-main sequence stars and to their accretion history. Outflows 
and jets are more powerful during the earlier stages and they 
tend to decrease in intensity as the star evolves toward a state 
of lower level of accretion. Unfortunately, there is little informa- 
tion in the literature about the stars of the sample showing neon 
emission consistent with outflows/jets. They are generally clas- 
sified as Class II sources, and the levels of disk mass accretion 
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rate are in the typical level of T Tauri stars ( ~ 10 8 M yr ';e.g., 
iGiidel et alJl2007t iGudel et~ai1l2010l and references therein). 

ISacco et aD (120121) also argue that the ratio between [Ne II] 
fluxes measured by Spitzer and VISIR are consistent for transi- 
tion disk and pre-transition disk sources. Our detected sample 
does not include targets with such a classification, but Class II 
stars (except V892 Tau) that show VISIR fluxes several times 
lower tha n Spitzer fluxes, co nsistently with the Class II stars ob- 
served bv lSacco et al .(2012). The latter argued that the emission 
in transition and pre-transition disk objects would be located 
within a few AU from the central star, in contrast with Class II 
sources that would emit [Ne II] from both the inner region and an 
extended envelope. In our three targets with low velocity shifts, 
extended emission is found for both FS Tau A and V892 Tau, 
whereas CoKu Tau 1 shows no evidence of extended emission 
(but its emission is either due to a bipolar j et or a disk) . Thus , 
overall, our data support the interpretation by Sac co et al. 
that extended emission is present in some Class II sources. 

In three stars the [Ne II] line could be interpreted with 
disk emission: CoK u Tau 1 (this work), AA Tau, and GM Aur 
dNajita et al.l 120091) . The first two are classical optically thick 
disks and GM Aur is a more evolved transitional disk. The first 
two also have high inclination angles: i = 87° ( CoKu Tau 1; 
Robit aille etaP l2007h and i = 75° (AA Tau. iBouvier et al.1 
1999b . whereas GM Aur has i = 54° (ISimon et all 12000). 
As discussed in this paper, the disk origin for CoKu Tau 1 
is difficult to as certain and could be attributed to the bipolar 
jet. Interestingly. ICox et al.l d2005l) r eports the presence of a jet 
in AA Tau. Although Naiita et al.1 (120091) interprets the broad 
(FWHM=70 km s -1 ) [Ne II] line as due to the disk, a similar 
situation as in CoKu Tau 1 could occur, and the emission could 
be due to a bipolar jet in AA Tau. However, in the case of GM 
Aur, its moderate inclination, its transitional disk status, and the 
absence of a reported jet all suggest that the emission may really 
come from the disk. In any case, the X-ray luminosities of the 
disk sources are of the order of 10 30 erg s _1 , sufficient to effec- 
tively heat the disk and to be responsible for the neon emission. 

Different studies suggest that soft X-rays can drive power- 
ful photoevaporative winds that can disperse the disk effectively. 
Observationally, these winds translate in emission lines show- 
ing a small blueshift. These winds are stronger, and are effec- 
tive once the mass accretion rate has fallen to values lower than 
10~ 8 M yr~' for an X- ray luminosi t y Lx ~ 2 x 10 30 erg s _1 . 
Out of 12 targets from Sacc o et all (120121) . ten show evidence 
of small blueshifts within 18 km s~'of the stellar velocity. Only 
two targets were classified as Class II sources, while the rest was 
classified as transition or pre-transition disk sources. In our sam- 
ple, FS Tau A and V892 Tau also show a blueshift, and three 
additio nal sources s how such blueshift s in iPascucci & Sterzikl 
(120091) . Interestinglv. lSacco et all (1201 2l) found evidence that the 
FWHM of the detected lines in their sample increased with in- 
creasing blueshift. Our sample does not show this trend, with 
the majority of lines with FWHM around 26-27 km s _1 , except 
for CoKu Tau 1 that is much broader and EC 92 which is nar- 
rower but slightly above the VISIR instrumental resolution of 
10 km s _1 . In fact, a plot of FWHM vs v pea k from both our and 
Sacco's data sets shows no evidence of a correlation. The high 
num ber of [Ne II] li nes consistent with photoevaporative winds 
bv lSacco et al.l (120121) could suggest that photoevaporation of the 
disk is an effective mechanism. However, the significant number 
of non-detections of [Ne II] in our study may also indicate that 
photoevaporative winds may not be as efficient as suggested in 
all young disk-bearing stars. In addition, our detected sources 
show different behaviors, from jets, photoevaporative winds, to 



disk emission. It may, nevertheless, be that photoevaporation is 
the most efficient mechanism for [Ne II] e mission in transitio n 
and pre-transition disk sources, as argued bv lSacco et al.l(l2012l) . 
Disk emission appears, on the other hand, rare and subject to 
uncertainties such as the disk inclination. 

5. Summary and Conclusions 

We presented observations of the [Ne II] line at 12.81 pm at 
high spectral and spatial resolution obtained with VISIR- VLT 
in a sample of selected pre-main sequence stars with previously 
detected neon emission with Spitzer. The emission line was de- 
tected and spectrally resolved in seven stars. We interpreted the 
profiles and shifts according to three emitting mechanisms: 

- Lines consistent with photoevaporation. These are the lines 
detected in FS Tau A and the Herbig Be star V892 Tau. 
They show small shifts with respect to the stellar velocity 
and FWHM -26 km s _1 . 

- Lines consistent with shocked material in outflows/jet. The 
line centroids are blueshifted at large velocities in four stars: 
MHO-1, SST042936+243555, V853 Oph, and EC 92. To our 
knowledge, there is no previous record in the literature of the 
detection of a jet in these stars. On the other hand, for sev- 
eral stars known to be jet driving sources we did not detect 
neon emission. It is possible that in these cases the orienta- 
tion and width of the slit used prevented us from detecting 
neon emission at large velocities. 

- Lines consistent with disk emission. This is possibly the case 
of CoKu Tau 1 since the line is centered at the stellar velocity, 
has a symmetric profile and large FWHM ~ 55 km s _1 . The 
line also presents wings extending toward velocities close to 
100 km s _1 . The emission of the fast rotating gas is likely 
coming from regions close to the star. However, an interpre- 
tation of the line due to the bipolar jet cannot be excluded 
and may, in fact, be better in view of the high [Ne II] lumi- 
nosity in CoKu Tau 1 . 

For three stars in the sample, optical observations with 
UVES allowed the determination of the radial velocity, a cru- 
cial measurement for determining the line shifts. The profiles of 
the forbidden lines detected in the optical spectra were compared 
with the profile of the [Ne II] line to better constrain the emis- 
sion mechanism. In general we found a good agreement between 
the profiles of the optical forbidden lines and the [Ne II] lines, 
indicating a common emitting mechanism. Studies that combine 
observations in the two bands need to be done in a more system- 
atic way. 

Combining our results with previous studies based on high- 
spectral and spatial resolution of neon, we attempt to extract 
more general conclusions on the emitting mechanism, keeping 
in mind that the sample is incomplete and biased toward large 
neon luminosities. For stars with neon emission compatible with 
disk emission and photoevaporative winds, a relatio n between 
Nell] and X-ray luminosities is difficult to see. ISacco et all 
20121) found a correlation between the [Ne II] line FWHM and 
the peak velocity, but our study d oes not confirm this trend. In 
fact a combination of ISacco et alJ 's and our data sets shows no 
trend. 

The mass accretion rates of the stars showing lines consistent 
with photoevaporation are in the range at which this mechanism 
clears out the disk. Many stars in t his group are classified as tran- 
sitional disks lPascucci & Sterzikl d2009h : ISacco et al.1 d2012l) . but 
two stars with neon emission compatible with photoevaporation 
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presented here do not belong to this class ( V892 Tau and FS 
Tau A). These two objects are close binaries with separations of 
5 and 36 AU for V892 Tau and FS Tau A respectively, and sur- 
rounded by circumbinary disks. The systems were not separated 
by VISIR observations. An interesting result from our study is 
the low detection rate of [Ne II] line s consistent with p hotoe- 
vaporation, in contrast to the study by Saccoetal 1 (120121) . Even 
more interestingly, CoKu Tau 1 shows evidence of broad line in 
[Ne II] and in other optical forbidden lines, but the origin of the 
emission is likely from its bipolar jet, although disk emission 
cannot be excluded, nor a contribution from a photoevaporative 
wind. 

A deeper knowledge of the stars studied is needed. As an ex- 
ample, for the sources with neon emission compatible with jet 
emission, very little information is published. Important param- 
eters such as mass accretion rate, or even the presence of out- 
flow/jet activity from other tracers is not available. Observations 
of the [Ne II] line can be used to study jets, photoevaporation, 
and disk emission. A larger sample is needed in order to iden- 
tify the preferred emitting mechanism for the different types 
of stars. Although ground-based observations of [Ne II] from 
young stars are challenging, our work demonstrates that they are 
crucial to constrain the kinematics and physical region of the 
emitting source, and, thus, to determine the emitting mechanism 
of the [Nell] line. 
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Table 5. Properties of the optical UVES lines for CoKu Tau 1, FS Tau A, and V853 Oph. The [O I] lines at /16300 and ,16363 of V853 Oph were 
fitted using the sum of two Gaussians to account for the low and high-velocity components (LVC and HVC), while flux and equivalent width were 
calculated for the sum of the two components. 



Line 


Wavelength 


Center 


FWHM 


Flux (xlO- 15 ) 




EW 




(A) 


(kms 4 ) 


(km s" 1 ) 


(erg s _1 cm~ 2 ) 




(A) 


CoKu Tau 1 


Ha 


6562.819 


10.2(0.8) 


103.5(0.8) 


5.86 (0.1) 




-40.0 


He I 


5015.678 






< 1.0 


> - 


-14.0 


He I 


5875.621 






< 0.01 


> 


-0.1 


He I 


6678.151 






< 0.01 


> 


-0.1 


Li I 


6707.760 


5.1(1.8) 


35.3(3.8) 


0.10(0.1) 




0.7 


[N II] 


6548.050 


1.3(0.8) 


85.4(3.9) 


0.41 (0.1) 




-2.8 


[N II] 


6583.450 


8.5(0.8) 


77.3(1.6) 


1.33 (0.1) 




-9.1 


[O I] 


5577.339 






<0.06 


> 


-0.1 


[0 I] 


6300.304 


4.6(0.8) 


58.2(0.9) 


3.21 (0.1) 




-21.9 


[O I] 


6363.777 


3.9(0.8) 


54.6(1.6) 


1.02(0.1) 




-7.2 


[S II] 


6716.440 


7.0(0.8) 


53.5(0.8) 


0.85 (0.1) 




-5.7 


[S II] 


6730.815 


7.6(0.8) 


57.4(0.6) 


2.01 (0.1) 




-13.7 


Fell 


5018.434 






< 0.1 


> 


-2.0 


Fell 


5169.030 






<0.05 


> 


-0.7 


Fell 


5197.577 






<0.07 


> 


-1.0 


Fell 


6432.680 


38.1(4.4) 


292.8(12.3) 


0.9(0.1) 




-5.3 


FS Tau A 


Ha 


6562.819 


-6.0(0.4) 


71.8(1.3) 


832(5) 




-16.1 


He I 


5015.678 






< 1.4 


> 


-0.2 


He I 


5875.621 


1.9(0.7) 


41.2(1.9) 


3.2(0.6) 




-0.5 


He I 


6678.151 






< 129 


> 


-2.5 


Li I 


6707.760 


4.9(0.2) 


30.8(0.4) 


20.7(5.2) 




0.4 


[N II] 


6548.050 


-1.2(0.2) 


30.6(0.6) 


47.2(5.2) 




-1.0 


[N II] 


6583.450 


-3.4(0.2) 


31.8(0.4) 


181(5) 




-3.5 


[O I] 


5577.339 


-0.6(0.8) 


29.1(2.0) 


3.8(0.6) 




-0.6 


[O I] 


6300.304 


-0.6(0.2) 


26.0(0.4) 


295(5) 




-5.7 


[O I] 


6363.777 


-1.3(0.2) 


25.1(0.5) 


101(5) 




-2.0 


[S II] 


6716.440 


-0.9(0.2) 


26.4(0.5) 


97.9(5.2) 




-1.7 


[S II] 


6730.815 


-0.0(0.1) 


24.6(0.3) 


181(5) 




-3.5 


Fell 


5018.434 






< 1.0 


> 


-0.2 


Fell 


5169.030 






< 0.2 


> - 


-0.03 


Fell 


5197.577 


18.2(0.6) 


25.2(1.7) 


3.2(0.6) 




-0.5 


Fell 


6432.680 






< 2.1 


> - 


-0.04 



V853 Oph 



Ha 


6562.819 


1.5(0.3) 


117.5(0.8) 


235.2(0.8) 


-30.3 


He I 


5015.678 


1.8(0.3) 


21.4(0.7) 


6.4(1.4) 


-0.5 


He I 


5875.621 


3.0(0.1) 


36.1(0.5) 


35.0(1.4) 


-2.6 


He I 


6678.151 


2.8(0.2) 


17.8(0.9) 


7.6(0.8) 


-1.0 


Li I 


6707.760 


5.6(0.2) 


22.0(1.6) 


3.1(0.8) 


0.4 


[N II] 


6548.050 






< 0.1 


> -0.02 


[N II] 


6583.450 


-41.7(0.8) 


34.3(1.9) 


< 1.7 


> -0.2 


[O I] 


5577.339 


-4.0(0.8) 


49.6(1.9) 


6.5(1.4) 


-0.5 


[O I] (HVC) 
[O I] (LVC) 


6300.304 


-32.1(0.1) 
-2.8(0.1) 


32.8(0.1) 
18.6(0.1) 


13.2(0.8) 


-1.7 


[O I] (HVC) 
[O I] (LVC) 


6363.777 


-24.0(0.1) 
-3.9(0.1) 


26.9(0.1) 
11.2(0.1) 


3.9(0.8) 


-0.5 


[S II] 


6716.440 






< 1.6 


> -0.02 


[S II] 


6730.815 


-37.1(0.5) 


26.0(1.1) 


1.9(0.8) 


-0.2 


Fell 


5018.434 


1.6(0.2) 


14.9(0.5) 


5.4(1.4) 


-0.4 


Fell 


5169.030 


0.9(0.1) 


13.4(0.2) 


6.5(1.4) 


-0.5 


Fell 


5197.577 






< 1.2 


> -0.1 


Fell 


6432.680 






< 0.1 


> -0.01 



Note: The error in the center was calculated considering the contribution of the error in the Gaussian fit and in the stellar radial velocity. 
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Appendix A: Additional Material 

A.1. VISIR transmission spectra 

Figure [AT]provides the VISIR transmission spectra for both our 
targets and the telluric standard. 

A.2. Radial velocities used in this study 

Since the stellar radial velocity is a key parameter in determin- 
ing the centroid of the emission lines, it is worth to explain the 
choice of radial velocities used in this study. 

For the stars located in the Taurus region with no mea- 
surement of the radial velocity available in the literature, we 
adopted th e average radial velocity of Taurus, v ra d = 16.0 + 
0.4 km s _1 (iBertout & Genovall2006l) . We also use this value for 
the Herbig Be star V892 Tau . The 12 CO 7 = 2-1 observations by 
IPanic & Hogerheiidd (120091) were likely dominat ed by extended 
emission from the cloud. A search through the iHerbig & Belli 
( 1988) catalogue within 1' of V892 Tau for radial velocities of 
good quality also argued in favor of a radial velocity around 
15- ^kms- 1 forV892Tau. 

In addition, the radial velocities used for two stars 
(L1551 IRS5 and EC 92) are reported in the literat ure using the 
frame of the Local Standard of Rest (LSR), vlsr (ICovev et al.l 
l2006h . In these cases the radial velocities were converted to v ra d 
using the expression v ra d = vlsr— (t/o cos I +V e sin /) cos b— Wq, 
where / and b are the coordinates of the star in the Galactic co- 
ordinate system, and U Q - 1 0.3 km s" 1 , V H = 15.3 km s _1 , and 
W Q = 7.7kms~ 1 come from lNiinuma et al.l (1201 ll) . 



barycentric motion of the earth at the moment of the observa- 
tions to obtain the final radial velocity measurement. The radial 
velocities obtained with this method are presented in Table IA.1I 
where we also included the radial velocities from the literature, 
when available. Errors are reported between parentheses. For 
one object (V853 Oph) there is no previous measurement of the 
radial velocity, while for FS Tau A there are two measurements 
available in the literature, both having uncertainties larger than 
5 km s _1 . For FS Tau A and V853 Oph the precision achieved 
in this work is < 1 km s _1 . Finally, due to the lower signal-to- 
noise ratio obtained in the spectrum of CoKu Tau 1, fewer pho- 
tospheric lines were detected. Therefore, the precision achieved 
in the radial velocit y is lower for this star (17.5 + 2.4 km s~') 
than that derived bv lWhite & Hillenbrand! d200l . although our 
value is consistent. Therefore we used their radial velocity in this 
paper. 

In addition, we obtained the vsin i for the three stars by 
comparing the obser ved spectra with rota t ionally broadened s yn- 
thetic spectra (e.g., iBertone et aD l2008t iMuller et aD 120 lib of 
r eff = 3850 K and logg = 4.5 for CoKu Tau 1 and FS Tau A. For 
V853 Oph the values used were T e ff = 3370 K and logg = 4.5, 
the lines are spectrally unresolved. The results are reported in 
Table |AT] 

A.3. Lines detected in UVES spectra 



A.2.1 . Determination of the stellar radial velocity from UVES 
spectra 

To determine the radial velocity of our program stars, we used 
the optical spectra obtained by UVES. Optical spectra are popu- 
lated by a series of absorption lines from the stellar photosphere, 
where the lines are shifted due to the stellar motion with respect 
to the observer, the radial velocity. In order to determine the line 
shift, each spect rum was compared w ith a high-resolution syn- 
thetic spectrum (IBertone et al.ll2008l) of a star having the same 
spectral type as the star observed. Where an appropriate syn- 
thetic spectrum was not available, we used a synthetic spec- 
trum of similar spectral type as the star observed. For this pur- 
pose we used t he in teractive IDL-based software described in 
ICarmona et al.l (l2010l) . This tool allows the synthetic and target 
spectra to be displayed, and to calculate the radial velocity of the 
observed spectrum using the cross-correlation technique. The 
maximization of the cross-correlation function is a widely-used 
procedure to determine rad ial velocities (e.g., iTonrv & Davisl 
Il979t lAllende Prietol l2007h . If we have two arrays, S corre- 
sponding to the stellar spectrum and T the template spectrum, 
the cross-correlation between the two will be a new array C de- 
fined by the following expression: 

Ci = Yj T k S k + i (A.l) 

A- 

The maximum value of the cross-correlation function C will 
correspond to the element i = p where p is the shift in pixels 
between both the stellar spectra S and the template spectrum T. 
The position of the maximum of the function (the radial veloc- 
ity) and its error were calculated with a Gaussian fit. The ra- 
dial velocity measured by cross-correlation was corrected by the 
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Table A.l. Radial velocities obtained from UVES spectra 



Name 


v sin i 


v radial 


Vradial 






(km s" 1 ) 


(kms- 1 ) 


(km s" 1 ) 






This work 


This work 


Literature 


Reference 


CoKu Tau 1 


10 (3) 


17.5 (2.4) 


15.0(0.8) 


White & Hillenbrand (2004) 


FS Tau A 


15(5) 


17.1 (0.1) 


24(7) 


Folha & Emerson (2000) 








18(5) 


Eisl6ffel&Mundt(1998) 


V853 Oph 


< 7 


-8.9 (0.1) 
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1 ' 6r MHO 2 (2011-01-21; 

L 

J*" 






1 - 8 CoKu Tau1 

1.6 




3.0 
2.5 
2.0 
1.5 
1.0 
0.5 



FS Tau A (2010-0 


1-04) 









1.2 
1.0 
0.8 
0.6 



L1551 IRS 5 




1.6 

1.4 

1.2 

1.0 

0.8 

0.6 

0.4 
0.2 



UY Aur A (201 1-01 -23) 




1.6 

1 1.4 

5 12 
| 1.0 

« 0.8 

o 0.6 

2 0.4 
0.2 



UY Aur B (201 1-01 -23) 







12.808 12.812 12.816 




12.808 12.812 12.816 
Wavelength (nm) 



1.6 
1.4 
1.2 
1.0 
0.8 
0.6 



EC 92 (2)09-05-31) 



12.808 12.812 12.816 



Fig. A.l. VISIR transmission spectra. In black the spectra of the star and in red the spectrum of the telluric standard. We included 
one spectrum for the stars observed in two days. The spectra are presented before the barycentric and radial velocity correction to 
visualize the telluric absorption lines and the observed science and standard star/asteroid spectra. 
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Fig. A.2. Lines detected in the UVES spectrum of CoKu Tau 1 . The spectrum was normalized to the continuum level. The fitted 
Gaussian profile is overplotted in green solid line. 
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Fig. A.3. Lines detected in the UVES spectrum of FS Tau A. The absorption feature next to [S II] (-167 16) is the absorption line of 
Cal 0*6717). 




Fig. A.4. Lines detected in the UVES spectrum of V853 Oph. The profile of [O I] /16300 was fitted with two Gaussians to account 
for the high and low-velocity components observed. 



